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C^ We propose an innovative scheme exploiting discrete diffraction in a two- 

^^ dimensional array of coupled waveguides to determine the phase and amplitude of 

CN the mutual correlation function between any pair of three telescopes of an astroint- 

' ^' erferometer. 
O 



The Van de Cittert theorem allows the determination of the angular shape of an incoher- 



^ ent source from the measurement of the complex spatial correlation function of the emitted 

c/2 light [1]. This feature is routinely used in astrointerferometry to achieve aperture synthesis 

>-C by the coherent combination of light collected from an array of telescopes aimed at the same 



astronomical target. In simple terms, high resolution images of stars are obtained from 



j> the Fourier transform of the measured complex correlation function (complex visibilities) 

O 

(^ which is sampled at frequencies determined by the projected baselines connecting pairs of 

O 

lO telescopes. The current challenge of astronomical interferometry is to enable simultaneous 



evaluation of the complex visibilities on a large number of baselines [2], to allow interfer- 
ometric imaging techniques of fast astronomical events such as exoplanet transits [3] or 
nova/supernova explosions [4]. 
rS To date, most of the existing or planned interferometric facilities are using bulk optics to 

perform multiple telescope beam combinations. While this solution helps to keep costs rela- 
tively low, the setups are typically complex, bulky and require special construction to keep 
mechanical vibrations to an acceptable level. Due to their inherent thermal and mechani- 
cal stability, integrated optic beam combiners [5] could reveal decisive to perform multiple 
beam combination effectively. Recently, a planar integrated chip allowing the measurement 
of the phase over 6 baselines among 4 telescopes has been realized and tested [6]. Due 
to geometrical constraints, the complexity of planar integrated optical components grows 



rapidly as more telescopes and baselines are added. Three dimensional geoemtries, such as 
those enabled by laser-writing techniques [7, 8], could relax geometrical constraints e.g. by 
allowing fiber cross-overs, resulting in a design simplification of multiple beam combiners. 

In this paper, we put forward an elegant and innovative scheme allowing the measurement 
of the spatial coherence function sampled by three telescopes by means of a three dimensional 
integrated component. The device exploits discrete diffraction effects in a two-dimensional 
array of coupled waveguides to perform a two-quadrature measurement of the phase delay 
between all possible baselines in a 3-telescope stellar interferometer (see scheme in Fig. 1(a)). 

Starting point of our work is the consideration that free space diffraction entails intrin- 
sically the capability of measuring the wavefront profile. Phase retrieval from Gabor-type 
holograms (also known as shadowgraphy) is probably the most striking example of exploita- 
tion of diffraction for phase metrology purposes [9, 10]. Light propagation in an array of 
evanescently coupled waveguides is similar to conventional diffraction, however it is bound 
to a system with finite degrees of freedom. 

Here we demonstrate that the measurement of the discrete diffraction intensity pattern 
at the end of a N=3x3 waveguide array can be uniquely related to the phase and amplitude 
of the mutual correlation function of three monochromatic fields coupled to suitable input 
waveguides. We chose to number of the sites of the 2D array of waveguides progressively 
from left-to-right and from top-to-bottom. Let Ai, A^ and Ay^ be the complex amplitude of 
the light of the three beams at the input face of the waveguide array and injected in the 
sites i, j and k, respectively. The field at the output of the n*^ waveguide of the array is a 
linear combination of the input fields: 

Un = an,iAi + OnjAj + an,k-4k. (l) 

The coefficients an,i, Onj and an,k are a function of z and can be calculated exactly [11] or 
obtained by integration of the field propagation equations over the sample length L: 

1 
where L^ is the coupling length of the fiber arrays [11] and Ck,i is the matrix describing the 
relative strength of the coupling between neighbouring waveguides. The resulting intensity 



pattern Jn =< |f/nP > can be cast in the following matrix form: 

N 



-^n = 2J '^n,k</k, (3) 



k=l 



provided that we define the Jk as: 

ji=<iAr>, j2=<iAr>, J3=<i^kr> 

J4 + ^ J5 =< AiA* >= Ty J4 - 2J5 =< ^jA* >= r*j (4) 

Je + i J7 =< AiAl >= Tik Je - iJr =< ^kA* >= T-k 
Js + iJg =< A^k >= Tjk Js - iJg =< A^^A* >= r*k 

Here the functions F = |r| exp(i$) are the complex field correlation functions of all possible 
pairs of incoming fields [1]. In particular, it is easy to recognize that the variables J4, Jq 
and Js represent the real part of the corresponding correlation functions, while J5, J7 and 
Jg are their imaginary parts. The exact retrieval of the mutual coherence of the input fields 
is possible whenever the matrix a„ ^ is invertible and well conditioned. We numerically 
calculated the matrices for all possible combinations of three input fields on a 3x3 square 
array of waveguides of varying length. The optimal input combination and array length 
allow obtaining a maximum error of 10% on the retrieved correlation function components 
{Jn} with a 1% precision measuerement of the output intensity fields {In}- Figure 2 shows 
a possible output of such a configuration for three fully coherent injected fields exploiting 
the indicated phase differences. 

Better insight in the potential of our interferometric scheme is given by the results of 
the simulated observation of a binary star. Figure 3(a) shows the amplitude of the spatial 
correlation function expected for the binary star depicted in the inset, and observed by the 
array of Figure 1. The traces of the sampled points during 6 hours of observation with the 
astrointerferometer are superimposed to the image. For simplicity, it has been assumed that 
the star is located along the direction of the rotation axis of the Earth, and the telescope array 
is located at the pole where the star is observable. For the chosen baselines, the amplitude of 
the visibilities and the closure phase (sum of the phases of the three correlation functions [12]) 
during the observation time are shown in Figure 3(b-c), along with the estimated values and 
measurement errors obtained from the realistic modeling of the discrete beam combiner. 
In the model, we calculate for the selected observation points the expected intensities at 



the output of the discrete beam combiner {In}- We then model the error in the intensity 
measurement by adding to each {/„} a pseudo-random number representing the photon shot 
noise of detected intensity. We assume that the detector has a minimum shot noise of 0.5%, 
which is typical for 16 bit detectors working near the saturation level, which we set equal 
to the maximum value of the {/«}• For each of the 100 noise realizations we calculate the 
retrieved values of {Jn}- These sets are then used to calculate the average value and the 
standard deviation of the amplitude of the visibilities and the closure phases. 

The simulations point out that the proposed beam combiner could reach a precision of the 
normalized visibility measurement of ±0.05, while the the closure phase could be determined 
with a error below ±0.5 degree. These values are in line with the current performance of 
existing interferometers, we point out that the performance of the proposed wave coherence 
meter can be improved, e.g. by the optimization of the geometry of the fiber array. 

In conclusion, we have demonstrated theoretically that a two dimensional array of coupled 
waveguides can be used to determine simultaneously the amplitude and the phase of the 
correlation function of a light source, sampled with an array of 3 apertures (telescopes). The 
device could be realized by means of laser writing techniques and can be virtually scaled 
up to an arbitrary number of telescopes and baselines. Applications to both fringe tracking 
and visibility measurement of astronomical targets are foreseen. Moreover, because of the 
small dimensions and lightness of integrate optics, beam combination by means of discrete 
optics could result attractive in future space-borne astronomical interferometers [13]. 
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List of Figure Captions 

Figure 1. A pictorial view of the concept of the integrated discrete beam combiner. The 
hght collected by three telescopes arrenged as in (a) is injected in three selected sites of a 
two dimensional array of coupled waveguides ((b) Left). The resulting discrete diffraction 
intensity pattern at the end of the array ((b) Right) gives unique information about the 
amplitude and phases of the mutual correlation function of the observed star along the 
three baselines of the telescope array. 

Figure 2. Example of an output resulting from the injection of three equally intense, fully 
coherent beams in the highlighted waveguides of a sample length of L = l.SL^ {L^. is the 
coupling length of the fiber array, for the definition see [11]). The phase differences between 
the various input fields are reported at the bottom of the picture. 

Figure 3. Simulation of the observation of a binary star, (a) Amplitude of the spatial 
correlation function of the light emitted by the binary star depicted in the inset. The traces 
of the projected baselines of the interferometric array during an observation of 6 hours 
are superposed to the plot. Parameters of the binay star: speraration = 7 mas, ratio of 
intensities = 2.512 (corresponding to 1 magnitude). Retrieved amplitude (b) and closure 
phase (c) of the correlation function sampled along the trajectories of the projected baselines, 
as simulated by our model of a noisy discrete coherence meter. The expected values of the 
correlation function are represented by the continous lines. 
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FIG. 1: A pictorial view of the concept of the integrated discrete beam combiner. The Ught collected 
by three telescopes arrenged as in (a) is injected in three selected sites of a two dimensional array 
of coupled waveguides ((b) Left). The resulting discrete diffraction intensity pattern at the end 
of the array ((b) Right) gives unique information about the amplitude and phases of the mutual 
correlation function of the observed star along the three baselines of the telescope array. 




FIG. 2: Example of an output resulting from the injection of three equally intense, fully coherent 
beams in the highlighted waveguides of a sample length of L = 1.8Lc (Lc is the coupling length of 
the fiber array, for the definition see [11]). The phase differences between the various input fields 
are reported at the bottom of the picture. 
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FIG. 3: Simulation of the observation of a binary star, (a) Amplitude of the spatial correlation 
function of the light emitted by the binary star depicted in the inset. The traces of the projected 
baselines of the interferometric array during an observation of 6 hours are superposed to the plot. 
Parameters of the binay star: speraration = 7 mas, ratio of intensities = 2.512 (corresponding to 
1 magnitude) . Retrieved amplitude (b) and closure phase (c) of the correlation function sampled 
along the trajectories of the projected baselines, as simulated by our model of a noisy discrete 
coherence meter. The expected values of the correlation function are represented by the continous 
lines. 



